Fluctuators in disordered metallic point contacts: a simulation approach 
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Electron transport through amorphous monatomic metallic structures generated earlier by molecu- 
lar dynamics simulations is studied numerically. The interference of electronic trajectories backscat- 
tered by the structural disorder probes the multistable structural relaxations responsible for low- 
frequency noise in real metallic contacts. The structure of these modes is visualized; the dependence 
of noise magnitude on size and structure of the modes is studied. The transition from multistable 
behaviour to a more complex one is observed for temperatures far below the melting temperature. 
The current fluctuations observed numerically resemble the complex behaviour reported earlier for 
current noise in small metallic structures at moderate temperatures. 



PACS 72.10-d, 72.70+m 

I. INTRODUCTION 

In recent years substantial progress in the study of the 
nicker noise problem was achieved. This is mainly con- 
nected with experiments on small conducting systems, 
where the 1// noise was decomposecffil into a sum of 
contributions due to "elementary" defects with internal 
degrees of freedom - "fluctuators" . These defects switch 
between two possible configurations and cause the "tele- 
graph" resistance noise in small contacts. In metallic 
systems they are believed to be related to structural de- 
fects which at low temperatures are seen as the well- 
known two-level systemsETEl. We would also like to men- 
tion works where the two-level fluctuation studied exists 
in atomic size ccjatacts produced with help of the break- 
junction methods. However, there is still no clear under- 
standing of the microscopic nature of these defects. 
Many experimental results can be explained!^ by the 
model of independent fluctuators described by rate equa- 
tions for the occupation numbers. The soft-potential 
approach which relates the fluctuators to some "soft" 
double-well effective—potentials, has been used to de- 
scribe this pictureQllj. It fails, however, to give a mi- 
croscopic insight into the problem. On the other hand, 
some experiments revealed a much more complex be- 
haviour exhibiting in particular fluctuator interactions 
and non-stationary chaotic behaviour of the fluctuator 
parameters! These features can not easily be described 
within the framework of the soft-potential model. Thus 
a microscopic analysis of the problem would be desirable. 
However, a description of disordered materials based on 
purely analytical methods meets obvious difficulties. Nu- 
merical simulations based on molecular dynamics (MD) 
methods are used successfully to model the structure of 
glassy materials as well as their thermodynamic proper- 
ties including the density of states. In particular, many 
details related to low-energy vibrational excitations in 
these materials have been revealed this wayEJ. To our 
knowledge, electron transport through disordered media 
has not been treated by these numerical methods. 
Real metallic point contacts studied experimentally can 



be as small as 3 — 10 nm in diameter and thus consist of 
only ~ 10 3 — 10 4 atoms which is just of the order accessi- 
ble to numerical simulations. Moreover, the structure of 
these constrictions is not expected to be perfectly crys- 
talline especially in view of the amorphous films used to 
prepare these point contacts. 

Therefore, it is attractive to use a simulation approach 
to study slowly relaxing metastable configurations in dis- 
ordered metallic structures with an application to low- 
frequency noise. In the following we will give the results 
of a first such study of noise within the framework of 
simplified model considerations. 

The main results of our paper described in the follow- 
ing are: 1) We have identified local atomic arrangements 
leading to metastable configurations. 2) A transition 
from a simple telegraph signal to a more complex be- 
haviour with increasing temperature was observed which 
reproduces the experimental behaviour; it is important 
that the ratio of this "critical" temperature to the melt- 
ing temperature agrees with experiment. 3) The magni- 
tude of the resistance fluctuations and its dependence on 
the number of atoms forming the "fluctuator" is consid- 
ered; the results of the simulations are in agreement with 
experiment. 



II. THEORETICAL ASPECTS 

Our aim is to analyze numerically charge transport 
through disordered— structures prepared previously by 
molecular dynamicaia The structures studied are consid- 
ered to model a 3D cubic microbridges (with side-length 
a) between bulk contacts. The atoms forming the struc- 
tures are taken to be scatterers for the incident electrons. 
We consider the weak scattering limit where the mean 
free path of the electrons, Z, is larger than a. This implies 
that the scattering cross-section s obeys the inequality 
riiS < a" 1 where rii = N/a 3 and N is the total number 
of atoms in the structure. In real glassy metals the mean 
free path is significantly larger than the electron wave- 
lengths, and therefore, the approximation is justified for 
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nanoscale point contacts. 

To study the conductance changes caused by "jumps" 
in the structural configurations we relate them to the 
changes of the interference contribution to the backscat- 
tered electron flux in the contact region. The weak scat- 
tering approximation allows us to neglect multiple scat- 
tering, restricting our calculations to a minimal number 
of scattering events. This corresponds to the concept of 
local interference which has been introduced, in Ref. |l3| for 
the case of scattering by complex defectsH. Mesoscopic 
effects in bulk samples caused by local interference of a 
more general type have been considered in Ref. [l5| 
The standard approach to the conductance G of small 
contacts exploits the Landauer formula (see e.g. Ref. [r]) 
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where R a is the reflection coefficient for channel (mode) 
a. In a simulation one needs to identify these modes a. 
This can easily be done for 2D quantum point contacts 
with smooth "adiabatic" boundaries!]] whereas for realis- 
tic structures one meets with considerable problems. In 
particular we mention Ref. [l^ where for the case of a 
short 2D channel the conductance is calculated numeri- 
cally by matching the "bulk" plane wave modes to the 
channel's transversal modes. 

Here, we will apply the "wave optics" approach which 
was developed for the transport in structures with a small 
number of scattering events and is in the spirit of Ref. 
|l8| . In Ref. [l9] this approach was developed and success- 
fully applied to analyze experimental data on mesoscopic 
transport in 3D ballistic point contacts with diffusive sur- 
roundings. We also want to mention the simulations p£ 
quasi-ballistic transport in static disordered structuresE2l 
based on the Green's function method whose results agree 
well with those of Ref. 

The "wave optics" approach exploits the fact that for 
metallic point contacts transport is always semiclassi- 
cal in the sense that the number of quantum channels 
is large, corresponding to the "geometrical optics" ap- 
proximation for the de Broglie waves. For an ideal point 
contact, separating the bulk leads, the modes a corre- 
spond to the "optical beams" cut by the orifice from 
the incident plane waves k. The "quantization" of the 
modes originates from the Fourier expansion of the ker- 
nel of the Fresnel-Kirchhoff integral which describes the 
diffraction of the incident plane wave by the aperture in 
an opaque screen, see e.g. Ref. The number of the 
modes, ~ (kpa) 2 (with kp the Fermi wave vector), is 
the number of independent solutions of the diffraction 
problem for a given |k| wkp. 

To calculate the backscattering coefficients R a within 
this approach we use the standard perturbation theory 
of scattering with Schrodinger's equations 2 ]. We consider 
the scattering of the incoming modes ipa which, as dis- 
cussed above, at distances of the order of the atomic scale 
can be approximated by plane waves with wave vectors 



k Q . The resulting wave, scattered by a scatterer i at 
point Hi, is 

r) = f a(R ; )F : • exp (ik\r- R, Q) . (2) 

Here Fi is the scattering amplitude which in the following 
we assume to be equal for all scatterers (F = (s/4-k) 1 ^ 2 ). 
Including multiple scattering would add higher order 
terms in F/a. For a random distribution of scatterers 
the sum over all possible configurations would give a fac- 
tor ~ TV 1 / 2 to the correction of the wave function Eq. 
|2|. Thus the total contribution of trajectories involving 
multiple scattering is by a factor s 1 / 2 A fl / 2 /(a(47r) 1 / 2 ) = 
s 1 / 2 ^ 2 'a 1 / 2 1 '{Ait) 1 / 2 « 1 smaller than the contribution 
due to single scattering, i.e. neglecting multiple scatter- 
ing is justified. 

The coefficient R a is by definition the ratio of the 
"backscattered" current to the "incident" current of the 
mode a. In our approximation we have: 
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Here ip s , a is the sum of the waves scattered by the scat- 
terers i, t>j_ iQ is the velocity of the incident electron 
wave perpendicular to the aperture. The apparent di- 
vergency due to vx,a in Eq. || is integrable due to 
v±. a oc {hp — ka^) 1 / 2 , with k Q _L the projection of k Q 
onto the orifice plane, and does not affect the final result 
of Eq. [|. S is the orifice area while r is the 2D radius- 
vector in the detector plane. Fig. |l| depicts schematically 
two backscattered trajectories which interfere in the de- 
tector plane. 

From Eqs. || and || we get for the change of the conduc- 
tance due to backscattering: 
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(cos(Z(OX, r - Rj)) + cos(Z(OX, r - R,-))) 

Here OX is the contact axis. Eq. |] is the basic input for 
the numerical studies described below. The advantage 
of the coordinate representation used in our approach is 
its direct contact to the simulation of electron transport 
through the structure. In our calculations we consider 
amorphous structures where all atoms act as scatter- 
ers. For relatively clean crystalline structures only defect 
atoms are responsible for backscattering. We are mainly 
interested in the changes of the backscattering due to 
changes of the structures and use all atomic positions as 
input for Eq. |4|. In the absence of defects one would ob- 
tain the conductance of purely ballistic contact from the 
Landauer formula of Eq. [I]. 
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III. COMPUTATIONAL DETAILS AND RESULTS 

With the help of Eq. |I] we carried out numerical studies of 
the backscattered current for amorphous and crystalline 
monatomic metallic structures and for amorphous sele- 
nium. We would like to note, however, that dealing with 
monatomic metallic glasses one meets difficulties due to 
their strong tendency to crystallize. This problem is 
known both from experiments on real metallic glasses and 
from computer simulations. One consequence is that the 
structural relaxations are mostly irreversible, preventing 
the formation of multistable structures supposed to give 
rise to the "fluctuators" . We will discuss these problems 
in more detail in relation to the results of our simula- 
tions. Additionally we have calculated the backscattered 
electron flux by using a more complex structure of amor- 
phous Se, which is much more stable against crystallisa- 
tion. The results of the MD-simulationE3€-3 have shown 
that both soft sphere and Se glasses behave similarly 
with respect to their dynamics which in both systems 
at low temperature exhibits local relaxations. Although 
electron transport in real Se-based structures is hardly 
similar to that in metals, we believe that the simula- 
tions for the corresponding model structures are in any 
case instructive. Indeed, one should keep in mind that 
"fluctuators" are thought to be related to defect regions 
where the structure might be even more complex than in 
a monatomic metallic glass. 

The metalti^^pwctures were generated in earlier MD- 
simulationsOl!3 ! Eil. The interatomic interaction in these 
systems is given by a soft sphere potential V(r) = 
e(a/r) 6 + A(r ja) 4 + B where the energy-scale is given 
in units of e and the length-scale in units of a. In the 
simulations we used the units e = a = m = 1, m being 
the atomic mass. The density is kept constant at pa 3 = 1 
and the nearest neighbour distance corresponds to about 
l.lcr. To reduce the numerical effort the potential is cut 
off at r c /er = 3.0 and shifted by a polynomial A(r/a) 4 +B 
with A = 2.54 ■ l(T 5 e and B = -3.43 • l(T 3 e. Both the 
potential and the force are zero at the cut-off. SmcIi a 
potential is appropriate to model metallic systemsEi 
As described by Laird and SchobenB to bypass crystal- 
lization the model glassy structures were produced from 
well equilibrated liquids by a rapid quench. With typical 
quench-rates of the order 0.25ItB/\/m<j 2 /e the configu- 
rations are cooled to roughly half the glass temperature 
T g and left there for more than 1200 molecular dynamics 
time-steps (MDS), i.e. about 60 typical vibration peri- 
ods, to stabilize the potential energy. The local minima 
of the structures were found by minimizing the poten- 
tial energy-, with a steepest descent/ conjugate gradient 
algorithmca. To have a reasonable statistics for the struc- 
tures consisting of 500 and 1024 atoms in all 60 and 21 
model glassy structures were produced, respectively. 
The structural changes of these metallic glasses were ob- 
served in a previous molecular-dynamics simulation!! 2 ]. 
To study the dynamics of these soft sphere glasses they 



were heated with help of the "Velocity- Ver let" -algorithm 
in stages to temperatures ranging from 0.05 to 0.125 T g . 
The molecular dynamics was done with periodic bound- 
ary conditions to reduce surface effects. After heating 
the glasses to the desired temperature, the temperature 
was kept constant by scaling the velocities of the atoms. 
At each temperature the dynamics of the systems was 
followed for 9000 time-steps, which corresponds roughly 
to 500 periods of a typical vibration. 
To observe more complex structures and to reduce fur- 
ther the "finite size" effects additionally 15 soft-sphere- 
glasses with TV = 5488 were generated from hot-melts and 
used to study relaxations in a greater detailnJ. These 
glasses were heated to temperatures ranging from 0.05 
up to 0.15 T g and left there for 90000 MDS (an order of 
magnitude longer than the smaller systems). 
During the molecular dynamics simulations at constant 
temperature the configurations move to new minima of 
the energy landscape due to local hopping processes over 
small energy-barriers. To detect new configurations on 
the energy hypersurface the displacements of the atoms 
from the metastable equilibrium positions IV n of atom n 
in configuration i are measured and the total displace- 
ment of the structure is defined as: 

Ai? 2 (i) = ^(R„(t)-Rj 1 ) 2 (5) 

n 

where R„ (t) is the actual position vector of atom n. If the 
total displacement of the atoms exceeds a cut-off value 
and the residence time of the atoms in the new posi- 
tions also exceeds a minimal period of at least three soft, 
low-frequency vibrations (several 100 MDS) the new po- 
sitions of the particles are accepted as a new minimum 
configuration. The cut-offs of displacement and resident 
time, respectively, are chosen such that spurious minima 
are avoided. This procedure is described in greater detail 
in Ref. H||. 

These minima are monitored during the molecular dy- 
namics runs. The corresponding relaxations from one 
minimum to another can be either reversible or irre- 
versible. The irreversible jumps are mainly due to the 
instability of the monatomic metallic glasses. This in- 
stability manifests itself in the fact that during the ob- 
servation time some of the samples partially crystallized. 
Such a crystallisation is connected with a global rear- 
rangement of nearly all atoms in the configuration. This 
effect depends strongly on the size of the samples used in 
the simulation. The larger the system, the more stable 
are the configurations. The largest systems show greater 
stability towards crystallization; just one configuration 
(out of 15 with N — 5488 atoms) crystallized during the 
observation time at a temperature T = 0.15T g . But the 
tendency to lower the potential energy prevails during 
the simulation. 

Therefore, the first jumps after heating up to different 
temperatures are related to non-local rearrangements of 
the glasses, i.e. jumps with large displacements involving 
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many atoms and leading to more stable configurations. 
These jumps are followed by relaxations located on only 
a few atoms. A measure for the localization of a jump is 
the effective mass 



(AR) 2 
IARLJ 
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where (AR) 2 = J2n(^n ~ R«) 2 * s the square of the total 
distance AR between two successive minimum configura- 
tions (called "initial" and "final" positions of the jump). 
R^* denotes the respective initial and final positions of 
atom n, |AR^ ax | is the maximal distance a single atom 
jumps in this relaxation. Since we take the masses of the 
atoms comprising the structures as equal, M c g counts the 
effective number of particles contributing to the hopping 
process. 

In an extended computer simulation it is found that the 
displacement of a single atom in the relaxation is only 
a fraction of the nearest neighbour distance (typically 
a tenth of a bond length), i.e. relaxations at low tem- 
peratures are local processes with small displacements 
of the atoms contributing to this thermally activated 
jump. However, even the small displacements of the 
single activated atoms exceed the vibration amplitudes. 
Since the typical total distance between two minimum 
configurations is about one or several nearest neighbour 
distancesE3, the effective mass is of the order of 10 to 100 
atoms contributing essentially to the relaxations. It is 
found that these entities of some ten particles involved 
in the hopping processes form chain-ilike structures which 
move collectively along these chainsca. 
Since we are interested in multistable configurations, we 
restrict ourselves to the relatively rare reversible changes 
of the configurations, which means that the initial con- 
figuration and the final one of successive jumps are iden- 
tical. A measure for the reversibility of successive jumps 
can be given by the quantity crri 
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where AR„ = (IV n — R^) is the displacement vector of 
atom n in the relaxation AR (R^ being, as above, the 
initial and the final positions of atom n in the relaxation 
AR), and AR and AR' are the total jump- lengths of two 

successive relaxations given by AR 2 = X^=i (^n)- F° r 
a reversible jump one has crr> — 1, i.e. the same atoms 
are contributing equally to both jumps. If the jumps are 
uncorrelated, crr? would be of the order of crri = N p /N, 
where N p is the number of atoms involved in the hopping. 
During the MD-runs of the 60 glasses with 500 atoms only 
two structures undergo reversible jumps, whereas for the 
glasses with 1024 soft spheres only one reversible relax- 
ation occurred. The effective masses of the reversible 
relaxations observed in soft-sphere-glasses with N = 500 
and 1024 atoms range from 6 to 16 atomic masses. At 
T re 0.1T g half of the glasses with N = 5488 settled 



during the observation time to configurations where re- 
versible jumps involving up to four local minima are ob- 
served. For the simulations with 5488 atoms reversible 
relaxations with effective masses from 9 to 50 atoms were 
found. 

In order to calculate the electron back-flow we use inci- 
dent electron waves with k- vectors of length k = n/a. 
This value determines the "discretization" of the numer- 
ical integration in Eq. ^. The summation over the modes 
in Eq. [| becomes a sum over discrete angles which we 
take as A<j) — 2ir/ka and A9 — n/ka, where a is the side- 
length of the amorphous sample, for the metallic glasses 
the side-length is a = iV 1 / 3 ^. Using different grids for 
the numerical discretization changes the results by a few 
percent only, i.e. they are not crucially sensitive to the 
specification of the modes. To escape numerical problems 
with vanishing denominators in Eq. ^ we integrated over 
r not in the front orifice plane but instead in a "detector" 
plane shifted by approximately one atomic distance from 
the orifice plane. The integration becomes a double-sum 
with x = iAx and y = jAy where Ax = Ay = 1/k. The 
detector plane is taken about 2 to 3 times larger than 
the orifice; the thus neglected backscattered current is 
less than 10 — 25%. As value of the scattering ampli- 
tude F, we take (typical for metals) kpF ~ 1. One can 
show that for a relatively small contact volume (small a 
and N) the weak scattering approximation (sN/a 2 < 1) 
holds with this assumption. 

To test our procedure we applied it to the simple case 
of a self interstitial atom in an fee metal. This so called 
dumbbell configuration is formed by replacing_pne atom 
by a pair of atoms aligned in (100) directions. The el- 
ementary jump process of this self interstitial atom con- 
figuration is to a nearest neighbour site and involves a 
change of orientation to a different cubic axis, i.e. the 
dumbbell is now aligned in (010) or (001) direction. Us- 
ing the soft sphere potential we construct configurations 
comprising 501 atoms with the dumbbell in the middle 
of the structure aligned in (100) and (001) directions. 
The coordinates of the atoms are then used as input for 
Eq. H to calculate the backscattering of these structures. 
For the configuration with the dumbbell in (100) direc- 
tion, i.e. perpendicular to the detector plane (see Fig. |l|) 
we find Wb = 0.536, and for the two other orientations 
(001) and equivalently the (010) alignment, parallel to 
the detector plane, we calculate Wb = 0.540. From this 
we deduce that the relative magnitude of the fluctuation 
due to the reorientation of the dumbbell is about 0.7%. 
In the following we present the relative magnitudes of 
the fluctuations caused by jumps and relaxations in 
monatomic metallic glasses and defect structures as de- 
scribed above. Since we are interested in the changes of 
the backscattered flow due to relaxations of the samples 
we determine the minima observed in the MD runs and 
calculate Wb for these structures corresponding to the lo- 
cal minima of the potential energy. We neglect the con- 
tribution of atomic vibrations to the electronic backscat- 
tering which will be similar for different configurations. 
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In Fig. H we show Wb for a metallic glass with N = 500 
atoms. After the initial rearrangement connected with 
a decrease in potential energy, the system was stable for 
more than 20000 MDS before a reversible jump occurred 
with a total jump- length of AR = 1.03a and an effec- 
tive mass M e ft — 16m. The relative magnitude of the 
backscattered current is about 0.6%. 
Fig. H shows Wb for a system with N — 1024 atoms. 
After approximately 6000 time-steps a reversible jump 
occurred with a total jump-length of AR = 0.44<r and 
an effective mass M e g = 6m. The fluctuation due to the 
structural relaxation is approximately 0.2% of the total 
interference contribution. 

To simulate more complex multistable atomic aggregates 
larger samples are needed. Indeed as we have discussed 
the smaller systems are relatively unstable thus prevent- 
ing extended multistable atomic aggregates. Therefore, 
we have also studied the current transport based on the 
results for a glass with N — 5488 atoms. To shorten 
simulations of the interference pattern we restricted the 
summations over the scattering atoms R in Eq. | to a 
subset N s « 700 of the atoms including and surrounding 
the active ones. 

The conductance change of this glass heated to 0.1T g 
and aged for 90000 MDS is shown in Fig. |. In Fig. |(a) 
we show the signal of the backscattering for the detec- 
tor placed perpendicular to the arbitrarily chosen x-axis 
of the structure, in Fig. [|(b) we measure the backscat- 
tered current in the plane perpendicular to the z-axis 
of the point contact, i.e. we have placed the detector 
parallel to the x-y-plane of the sample of the configu- 
ration and the bulk leads are now perpendicular to the 
x-y-plane of the nanoscale point contact. The results 
of the backscattering are much more complex than the 
ones described above for the smaller systems: the com- 
plex behaviour of the conductance change shown in Fig. |] 
results from 15 jumps between eight minima denoted by 
A,B,C,...,H according to their occurrence during the sim- 
ulation. After an initial drop in energy the glass relaxes 
to a (meta)stable region in the configurational space, 
where we monitored these strongly correlated jumps be- 
tween 8 minima of nearly equal energy. After about 
17000 MDS a reversible jump (E^F^E) occurred with 
a total jump-length AR — 0.72c and an effective mass 
M cS = 16m. At t w 42000 MDS the glass returns to 
minimum A and relaxes to minimum C (AR — 1.22a 
and M e g = 18m) and again back to minimum A. After 
visiting two other minima G and H, the jump pattern 
(A^C^A) is repeated. 

The magnitude of the relative backscattered electron flux 
is about 0.037% for the contribution of the first reversible 
relaxation (E^F— > E) measured in a detector plane per- 
pendicular to the x-axis of the sample (see the experimen- 
tal setup in Fig. [j]) and 0.032% for the second pronounced 
fluctuation (A^C^A). 

The signal of these fluctuators changes drastically when 
we measure the backscattered current in the plane per- 
pendicular to the z-axis of the configuration. As de- 



scribed above we have placed the detector parallel to the 
x-y-plane of the sample. The change of the experimental 
setup enables us to detect inhomogeneities in the struc- 
ture and in the dynamics of the nanoscale point contact. 
The first reversible jump (E^F— ► E) causes a change of 
the relative magnitude of the backscattered electron flux 
of approximately 0.11% and the second one (A^C— >A) 
leads to a relative change in the current of about 0.66%. 
These changes of the signals with respect to the axis may 
be due to the strong displacements of the atoms in z- 
direction in both reversible relaxations. 
Additionally we calculated the backscattered current also 
for a partially crystallized defect structure of N = 500 
atoms. After heating to 0.05T ff the configuration begins 
to crystallize. During the observation time the sam- 
ple relaxes to another defect configuration. The ini- 
tial relaxation is accompanied by a large drop in energy 
(AE = 30.28e) while the following jumps lead to min- 
ima with similar energies (AE < 10~ 7 e), in the last re- 
laxation the system jumps to a minimum with slightly 
higher energy than the one preceding (AE = 0.00670e). 
The corresponding jump-length is AR — 13.70(7 for 
the initial relaxation. For the three following jumps 
(AR = 5.23(7, 7.75cr, 6.83cr) with small energy changes 
we find effective masses M c s — 68m, 74m, 72m, respec- 
tively. For the last jump the effective mass is about 74m 
and the jump-length is about 6.53cr. The jumps follow- 
ing the initial decrease in energy are strongly correlated, 
i.e. mostly the same atoms are involved. Fig. || shows 
the results for the calculations of the backscattering flow 
Wb plotted vs. the observation time. The amplitude of 
the fluctuations after the initial step is about 3%. 
So far we have considered monatomic metallic glasses, 
where due to crystallization effects long observation times 
are not accessible. Compared to the soft sphere glasses 
the models for selenium are more stable, and longer ob- 
servation times are possible without an appreciable drop 
in potential energy. Therefore, we also report results of 
the backscattered electron flux, normalized with respects 
to the incident flow, for a selenium glass. Although se- 
lenium is not a metal, the complex structure comprising 
chains and rings may be a good model for atomic entities 
forming the "fluctuators" . 

The interaction of the Se-aioms is described by a short- 
range three-body-potentialEj. The glasses are produced 
by quenching corresponding liquids with quench-tates of 
f = 10 12 - 10 14 K/s and pressures up to 10 GPaH Note 
that if the quenching would be too rapid density fluc- 
tuations of the liquid could be frozen in. One possible 
way to avoid the formation of resulting micropores is to 
apply lower quench-rates and/or pressure. The glasses 
are heated to temperatures of about 0.30T g and then 
observed for approximately 0.5 ns at every temperature 
(this corresponds to roughly 3400 high frequency vibra- 
tions). In Fig. ^ we plot the backscattered flow Wb ver- 
sus time for a glass of N = 750 atoms at T « 0.04T g . 
One can clearly observe a telegraph signal caused by a 
reversible jump occurring four times during the observa- 
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tion period. Its effective mass is M e s = 10 and its total 
jump-length AR = 0.5A. The amplitude of the fluctua- 
tion is about 0.42% of the total interference contribution 
to the backflow. 



IV. DISCUSSION AND CONCLUSIONS 

We calculate the conductance fluctuations via the stan- 
dard Landauer formula from the change of backscattering 
by the atoms comprising the junction. For the small junc- 
tions considered here electron scattering can be treated 
ballistically in the weak scattering limit. 
In geometries where multiple scattering is important dif- 
ferent methods have to be used. Todorov and Sutton, and 
Bratkovsky et al. combine a molecular dynamics simula- 
tion with a T-matrix approach in tight binding approx- 
imation to calculate the conductance .through atomic- 
scale metallic contacts during pullofEJC3. In agreement 
with experiment they observe conductance jumps due to 
atomic rearrangements in the tip. 

We want to point out that our approach is not suitable 
to reproduce the limiting case where the point contact 
merges into the crystalline structure. In this limit the 
conductance is no longer dominated by the scattering of 
the atoms in the junctions and the aperture. 
Let us now consider the characteristic magnitude of the 
non-stationary current fluctuations. As seen from Eq. ||, 
one can discriminate between two terms: (i) "classical" 
backscattering related to the diagonal terms in the sum 
of Eq. ^ which is independent of the coordinates of the 
scatterers and thus not sensitive to relaxations; and (ii) 
the interference contribution due to the N(N — 1) non- 
diagonal terms proportional to cosines and depending on 
the positions of the atoms. For a given mode k and for 
a given position r in the detector plane the standard de- 
viation of the interference term is ~ \J N(N — 1) ~ N. 
However, additional averaging over the - N 2 / 3 differ- 
ent modes and over the different positions of r reduces 
the standard deviation to N/(N 2 / 3 ) 1 / 2 (N 2 / 3 ) 1 ^ 2 = N 1 / 3 . 
Here positions separated by distances larger than ~ 1/k 
are to be considered as independent, concerning the in- 
terference, giving in our units about N 2 / 3 independent 
positions. Therefore, the total standard deviation of 
the interference contribution due to a displacement of 
all scatterers is smaller than the average backscattered 
flow by about a factor of N/N 1 / 3 — TV 2 / 3 . If one con- 
siders the current transport according to the concept of 
"quantum channels" , where in the weak scattering limit 
each elementary cell in the point contact orifice roughly 
corresponds to one channel, the standard deviation of 
the interference contribution is of the order of the contri- 
bution of one channel, in agreement with the picture of 
Universal Conductance Fluctuations (see e.g. Ref. |34| ). 
For a shift of a single scatterer (affecting N terms in the 
summation over i,j in Eq. ^) the effect is expected to 
be reduced by N 1 / 2 while for an uncorrelated jump of 



M atoms the expected effect is ~ M 1 ' 2 /N 2 ' 3 ^/ 2 ~ 
M 1 / 2 N~ 7 / & (normalized with respect to the average 
backflow current). Thus a jump of about 15 atoms out 
of 500 should give a fluctuation of relative magnitude 
- 0.003. 

As seen, the order of magnitude of the experimental re- 
sults is in agreement with this estimate. Some larger 
effects shown in Fig. || can be related to the fact that 
actually the jumps cannot be considered as incoherent 
shifts of M atoms; the coherency of the atomic motion 
can increase the observed effect. Qualitatively, the re- 
sults show the decrease of the relative fluctuation with 
increasing N and with decreasing M. The magnitude of 
the fluctuations gives information about the size of the 
relaxing region and about the character of the motion. 
The occurrence rate of the fluctuations is many orders 
of magnitude smaller than characteristic vibrational fre- 
quencies of the structure and thus we are obviously deal- 
ing with slowly relaxing degrees of freedom. However, 
for the simulated structures the rate is still much larger 
than a typical frequency for the telegraph noise observed 
in real experiments™]. The reason is again the instability 
of the structures studied with respect to crystallization 
which simply prevents long observation times which is 
in particular related to the small sizes of the systems 
and to the role of the periodic boundary conditions. In 
real systems relaxations would be slowed down due to 
strains exerted by the embedding material. In our point 
of view the simulations give a good model for the slowly 
relaxing multistable entities in real structures. The struc- 
tures in real point contacts will be somewhere in between 
the amorphous and partially crystalline state. In this re- 
spect we would especially like to emphasize that the mul- 
tistable defects were also demonstrated for the simulated 
crystalline defect structure of Fig. |5[ 
The simulated current fluctuations exhibit a rather com- 
plex time behaviour. The complex behaviour shown in 
Fig. ^ results from jumps connecting eight different con- 
figurations of nearly equal energy. Since the initial and 
final configurations are identical all jumps are reversible. 
Some jumps might trigger each other, others will occur 
randomly. E.g. after 42000 MDS a jump followed im- 
mediately by the return jumps occurs, 15000 MDS later 
the same jump is followed by two different jumps before 
returning to the start configuration. Actually the simu- 
lated picture resembles the one observed for real metal- 
lic contacts 2 ! at high enough temperatures referred to as 
"fluctuator melting". We would like to emphasize that 
the slowing down of relaxations with temperature pre- 
vented us from studies of the extremely low temperatures 
where simpler two state telegraph signals are expected. 
We want to emphasize, however, that the temperatures of 
our simulations are still much below the glass-transition 
and melting-|temperatures. We see a similarity to real 
experimentsa where "fluctuator melting" was observed 
at temperatures of about 100 K, i.e. much below the 
melting temperature of the metal. 

Our simulations show that typically the "fluctuators" in- 
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volve large groups of atoms comprising up to 100 parti- 
cles arranged mainly along some preferential direction. In 
Fig. fj] we show such a fluctuator identified by the of MD- 
simulations for the current noise shown in Fig. ^. The 
atoms contributing strongly to the relaxation are shown 
with their respective initial and final positions marked as 
full and shaded circles. The structure formed by the par- 
ticles resembles a chain with side branches. Note that the 
displacement of the atoms from the initial minimum to 
the second one is nearly one-directional (along the chain). 
This fluctuator corresponds in Fig. |^ to the pronounced 
reversible relaxation after 42000 MDS with an effective 
mass M e ff = 18m and a total displacement AR = 1.22er. 
Multistable configurations are more easily formed in the 
more complex Se structure. To some extent this hints 
towards a possible role of complex defects in the forma- 
tion of fluctuators. In Fig. |§| we show the atoms, which 
strongly participate in the telegraph noise shown in Fig. 
H at t < 250 ps. The complexity of the relaxation is 
reflected by the fact that atoms of two regions of the sys- 
tem strongly participate in the jump. One of the struc- 
tures of the relaxing entity is nearly one-dimensional and 
the atoms comprising this structure move cooperatively 
along this preferential direction. In particular, one can 
speculate on the apparent similarity between the helical 
structure in Se and the structure near the core of a screw 
dislocation. In this respect we would like to mention the 
experimental studiesu where the observed 1// noise was 
shown to originate from atomic motion near dislocation 
cores. 

One notes that the picture of the current behaviour is 
more complex than simple fluctuations which average out 
over long enough time scales which are not accessible in 
MD-simulations. We observe a number of relaxations 
and related current fluctuations with different relaxation 
times. The time window of our simulation is too narrow 
to determine whether these fluctuations obey a 1/ f law 
in the large ensemble or long time limits. 
To conclude, in the framework of a simple model we have 
done a first simulation of multistable current noise in dis- 
ordered contacts. The slowly relaxing modes responsible 
for the noise are visualized. The dependence of the noise 
magnitude on the size and structure of the modes was 
studied. In particular, the role of coherent atomic mo- 
tion was emphasized. An important feature is a dras- 
tic change in the fluctuator dynamics (a transition to a 
complex multistable character of fluctuations) at temper- 
atures of about a tenth of the melting temperature of the 
sample; this behaviour agrees with "fluciuators melting" 
reported for real metallic point contacts! 
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FIG. 1. Schematic plot of backscattered trajectories each 
involving one atom acting as scatterer of the incident electron 
wave. The incident wave is represented by its wave vector 
k. Filled circles depict some of the atoms in the point con- 
tact. Ri and Rj denote the positions of the atoms and r the 
two-dimensional radius vector in the detector plane. 



FIG. 5. Fluctuator behaviour of a configuration with 
N = 500 atoms, which partially crystallized at a tempera- 
ture T = 0.05T 9 . The last three jumps with effective masses 
M e H = 74m, 72m, 74m have jump-lengths A_R = 7.75a, 6.83a 
and 6.53a, respectively. These relaxations are strongly corre- 
lated with Cb.ri = 0.958, 0.996. Wb is plotted against time in 
units of MDS. 



FIG. 6. Fluctuator behaviour of a configuration with 
N — 750 Se-atoms, at a temperature T = 0.04T 9 . During 
the time-range t = 15 — 260 ps the glass switches between two 
configurations. The effective mass of the reversible jumps is 
M e fl ~ 10 and the jump-length A_R « 0.5A. The signal of the 
fluctuation is about 0.42% of the total contribution. Wt, is 
plotted against time in units of ps. 



FIG. 7. Reversible relaxation at T 
glass with TV = 5488 atoms. Full spheres depict the initial po- 
sitions of the atoms, shaded spheres are the atomic positions 
after the relaxation. Shown are all atoms with a displacement 
more than 0.4 of the maximal contribution to this relaxation. 
The effective mass of this relaxation is M e g — 18m and the 
total jump-length is AR = 1.22a. The largest contribution of 
a single atom to this jump is A7? max = 0.29a. The backscat- 
tering contribution of this jump is shown in Fig. ^| 



FIG. 8. Reversible relaxation at T 
a Se-glass with TV = 750 atoms. Full spheres depict the ini- 
tial positions of the atoms, shaded spheres are the atomic 
positions after the relaxation. Shown are all atoms with a 
displacement more than 0.2 of the maximal contribution to 
this relaxation. The maximal contribution to this relaxation 
is Ai? max = 0.44 A. The backscattering contribution of this 
jump is shown in Fig. ^. 



FIG. 2. Fluctuator behaviour of a configuration with 
TV = 500 atoms. Wb is plotted against time in units of MDS. 
The reversible relaxation at t m 25000MDS has an effective 
mass M e ff = 16m and a total jump-length AR = 1.03a. 



FIG. 3. Fluctuator behaviour of a configuration with 
TV = 1024 atoms. Wb is plotted against time in units of MDS. 
The reversible relaxation at t ~ 6000MDS has an effective 
mass M e ff = 6m and a total jump-length AR = 0.44a. 



FIG. 4. Fluctuator behaviour in a glass of 5488 soft spheres 
undergoing fifteen jumps between eight minima denoted by 
A,B,C,...,H. Different directions of detector planes are shown: 
(a) plane _L x-axis and (b) plane _L z-axis. Two direct re- 
versible relaxations (with effective masses M e s = 16m and 
M e g = 18m and corresponding jump-lengths AR = 0.72a 
and AR — 1.22a) can be observed. Wb is plotted against 
time in units of MDS. 
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